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We present new (2004 July) G750L and G140L Space Telescope Imaging 
Spectrograph (STIS) data of the Ha and Lya emission from supernova remnant 
(SNR) 1987A. With the aid of earlier data, from Oct 1997 to Oct 2002, we track 
the local evolution of Lya emission and both the local and global evolution of Ha 
emission. We find that the average Ha intensity has increased locally by a factor 
~ 3 for both blueshifted and redshifted emission, over periods of about 5 and 
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over about 5 years, for both components. The most recent observations allow us, 
for the first time, to directly compare the Ha and Lya emission from the same 
slit position and at the same epoch. Consequently, we find clear evidence that, 
unlike Ha, Lya is reflected from the debris by resonant scattering. In addition 
to emission which we can clearly attribute to the surface of the reverse shock, 
we also measure comparable emission, in both Ha and Lya, which appears to 
emerge from supernova debris interior to the surface. New observations taken 
through slits positioned slightly eastward and westward of a central slit show a 
departure from cylindrical symmetry in the Ha surface emission. No obvious 
asymmetry is seen in the interior emission. Using a combination of old and new 
observations, we construct a light curve of the total Ha flux, F, from the reverse 
shock, which has increased by a factor ~ 4 over about 8 years. However, due to 
large systematic uncertainties, we are unable to discern between the two limiting 
behaviours of the flux — F oc t (self-similar expansion) and F oc t 5 (halting 
of the reverse shock). Such a determination is important for constraining the 
rate of hydrogen atoms crossing the shock, which is relevant to the question of 
whether the reverse shock emission will vanish in < 7 years (Smith et al. 2005). 
Future deep, low- or moderate-resolution spectra are essential for accomplishing 
this task. 

Subject headings: circumstellar matter — shock waves — supernovae: individual 
(SN1987A) — supernova remnants 

1. INTRODUCTION 

SNR 1987A serves as a unique astrophysical laboratory for the study of the physics of 
shocks. Prior to the supernova (SN) explosion, the progenitor, Sanduleak — 69°202 (Type 
B3 I blue supergiant; Rousseau et al. 1978), was surrounded by a circumstellar ring system, 
consisting of an inner equatorial ring and two outer rings (Cassatella et al. 1987; Fransson et 
al. 1987; Crotts, Kunkel & McCarthy 1989; Crotts & Kunkel 1991; Wang & Wampler 1992; 
Crotts, Kunkel & Heathcote 1995). The impact of the ejecta with the circumstellar rings 
has been anticipated since the discovery of the circumstellar gas and various predictions had 
been made regarding the time of impact (Luo & McCray 1991; Luo, McCray & Slavin 1994; 
Chevalier & Dwarkadas 1995). As the blast wave propagates through the ambient medium, a 
double shock structure is established (Chevalier 1982). A shocked H II region resides behind 
the blast wave, and a contact discontinuity separates this region from the shocked ejecta. Ha 
and Lya emission from the reverse shock result from electron and ion impact excitation of 
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neutral hydrogen atoms that cross the reverse shock and enter the shocked ejecta (Borkowski, 
Blondin & McCray 1997, hereafter B97). As the hydrogen atoms are freely streaming, there 
is a unique mapping of the Doppler shift, Au, of the emitted photons to the depth of the 
supernova debris: 

Vq Ct 

where u is the initial frequency of the photon, the depth, z, is measured from the midplane 
of the debris along the line of sight, and t is the time since the supernova explosion. Since 
the equatorial ring is at an inclination of ~ 42° — 44° (Panagia et al. 1991; Burrows et al. 
1995; Plait et al. 1995; Sugerman et al. 2002), the blueshifted and redshifted emission are 
associated with the northern and southern sides of the debris, respectively. The emission 
mechanism is similar to the one present in Tycho's remnant (Chevalier & Raymond 1978). 

Since the work of Sonneborn et al. (1998, hereafter S98), the SAINTS (Supernova 
1987A INTensive Study) group has been observing the reverse shock emission with the 
Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). 
S98 detected the first definitive sign of impact between the blast wave and the equatorial 
ring, in 1997 April STIS spectral images in Lya. Michael et al. (1998a) built a simple 
kinematic model, based the data obtained by S98, and found that the high-velocity Lya 
emission could be explained by ~ 15,000 km s -1 hydrogen atoms crossing a reverse shock 
in the shape of a slightly prolate ellipsoid. Subsequently, Michael et al. (1998b) found that 
a model in which the reverse shock emission was confined primarily to an equatorial band, 
as proposed by B97, agreed better with new (1997 Sept & Oct) data. They also detected 
the presence of residual Lya emission that appeared to come from a volume interior to the 
reverse shock surface, and a departure from axisymmetry in both Ha and Lya emission, 
which correlated with observations of non-thermal radio emission (Gaensler et al. 1997). 
Michael (2000) and Michael et al. (2000) followed up with more STIS observations and 
proposed that "hot spots" are formed where the blast wave overtakes density protrusions 
on the equatorial ring. Pun et al. (2002) analyzed and modeled the optical and ultraviolet 
spectrum of the first hot spot, also known as "Spot 1" and found that the emission from 
the shocked gas came from slower (< 135 km s _1 ), radiative shocks. Lawrence et al. (2000) 
and Sugerman et al. (2002) analyzed the emergence of several more hot spots. Michael et 
al. (2003, hereafter M03), analyzed high- velocity (~ ±12, 000 km s _1 ) Ha and Lya emission 
profiles to map the geometry and development of the reverse-shock surface and found its 
shape to be consistent with a model in which the supernova debris expanded into a bipolar 
nebula. M03 also found evidence for the resonant scattering of Lya photons within the 
supernova debris and detected emission in both Ha and Lya that appeared to come from 
inside the reverse shock surface. 
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With the demise of STIS, it is an appropriate time to consolidate what we have learnt 
about the reverse shock emission and its evolution. In §2, we list the observations and 
describe our data reduction procedures. In §3, we analyze the evolution of the Ha and 
Lya emission from the reverse shock. Because telescope time was limited, we were unable 
to completely map the emission at each epoch. Therefore, the observations are sparse and 
we can directly compare observations at different epochs only in restricted regions where 
the observations overlap. In our most recent (July 2004) set of observations with STIS, we 
obtained spectra of the reverse shock in both Ha and Lya with the same slit. We compare 
these spectra in §4 and §5. In §6, we construct the light curve of the reverse shock Ha 
emission from new (July 2004) and old (Oct 1997 to Oct 2002) STIS observations. We 
also include ground-based data from Smith et al. (2005, hereafter S05), obtained during 
the commissioning run of the Low Dispersion Survey Spectrograph-3 (LDSS3) on the Clay 
Telescope of the Magellan Observatory. In §7 and §8, we interpret and summarize our results. 



2. OBSERVATIONS & DATA REDUCTION 

We analyze STIS observations of SNR 1987A, taken with the G750L and G140L gratings, 
from 1997 Oct 6 to 2004 July 23 (days 3878 to 6360 since the SN explosion). In Table 1, 
we list these observations ("pX/vY" means "Proposal X, Visit Y") and in Fig. 1, we display 
the corresponding slit positions overlaid on images of the circumstellar ring. All velocity 
scales are relative to the systematic velocity of SNR 1987A, v^ja = 289 km s _1 (Crotts & 
Heathcote 2000). 



2.1. CLEANING 

We employ standard techniques to clean the data, as described by M03. Each obser- 
vation consists of three or four exposures taken at dithered positions along the slit. The 
spectral resolutions are ~ 450 km s _1 and ~ 300 km s^ 1 for the G750L and G140L data, 
respectively. We reduce the two-dimensional spectra using standard calibration files. We 
remove hot pixels by utilizing the cr_reject procedure first written by Hill et al. (1997), using 
a user-supplied error cube and default settings for the clipping sigmas. We then combine 
the exposures. The resulting images (Fig. 2) contain not only spectra of the reverse shock 
but also a number of contaminants as indicated. 

In the G750L spectra, we see emission lines from the circumstellar ring at Ha, [N il] 
AA6548, 6584, [O i] AA6300, 6364, [S il] AA6717, 6731, and He I A6678. These are due 
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to residual recombination and cooling following the initial ionizing flash of the supernova 
(Fransson & Lundqvist 1989). We also see medium- velocity (~ ±3000 km s _1 ) emission 
from the radioactive core (defined as the central ejecta powered by the decay of radioactive 
species) of the supernova in Ha and [O i] AA6300, 6364. We mask out these contaminants and 
remove both background emission and stellar continuua using average profiles. We deredden 
the spectra using the extinction curve towards Star 2 (Scuderi et al. 1996). 

In the G140L spectra the main contaminants are N v AA1239, 1243 emission, broad 
damping wings from interstellar H I absorption, geocoronal emission at Lya and a few 
fainter lines, and interstellar line absorption due to C I AA1191, 1259 and N I A 1199. We 
remove the geocoronal emission by subtracting an average spectrum and then masking out 
the central slit. We mask out the N v AA1239, 1243 emission, since its intrinsic distribution 
from the ring system is unknown. We construct composite absorption profiles by spatially 
integrating and combining spectra from many separate observations and use these to correct 
for the interstellar absorption lines. To correct for broad damping wings due to interstellar 
Lya absorption, we divide the data by an absorption profile corresponding to Nhi = 3 x 10 21 
cm -2 , following M03. In all of the observations, the statistical errors are negligible (~ 0.1%). 

2.2. COMPONENTS OF THE REVERSE SHOCK EMISSION 

The Ha emission streaks labeled RS in the G750L observation in Fig. 2a are evidence 
that the slit is intersecting a curved surface of the reverse shock, on which the emitting 
hydrogen atoms have a unique velocity at every position. We call these features "surface 
emission". In addition to these streaks, one can also see faint Ha emission extending to 
lower velocity at the same vertical positions. If this emission comes from freely-streaming 
hydrogen atoms in the supernova debris, then it must originate from within the reverse 
shock surface. Therefore, we call it "interior emission" . When we refer to the "reverse shock 
emission", we mean both the surface and the interior emission (see §5 and §7.2). Wherever 
possible, we isolate the former from the latter. We emphasize that the interior emission is 
distinct from the core emission (see §2.1). In the G140L Lya spectra (e.g., Fig. 2b), the 
interior emission is so prominent that it is difficult to distinguish from the surface emission. 
As noted by M03, the interior emission in Ha is not subject to resonant scattering, as is the 
case for Lya. Therefore, we use the Ha surface emission streaks as templates to locate the 
surface emission in the Lya spectrum (see §4). 



- 6- 



2.3. SPATIAL-SPECTRAL CORRECTION 

While the emission from the SN core fades, that from the reverse shock brightens. Since 
the G750L 2" observations were taken at earlier epochs, reducing them is more challenging 
than for the later observations, due to the low signal-to-noise ratio and the presence of 
multiple contaminants. Furthermore, in STIS observations taken with the 2" slit, each point 
along a vertical (i.e., spatial) cut corresponds to a different velocity, owing to the curved 
shape of the reverse shock surface viewed through the slit. To correct for this curvature, we 
apply a spatial-spectral correction to morph each feature into a "bar" , using the [O i] A6300 
line from the unshocked ring as a reference (Fig. 3). This procedure not only facilitates the 
masking out of contaminants, but also ensures that the spatially summed data are placed in 
the correct velocity bins. The correction is applied for the (X'5 observations as well. 



3. LOCAL EVOLUTION OF REVERSE SHOCK EMISSION 

Since no two subsequent observations with the (X'2 slit have the same slit positions, we 
can only measure the evolution of intensity in the areas of overlap, which are in the form of 
parallelograms, labeled HI to H6 (Fig. 4) and LI to L6 (Fig. 6), for G750L and G140L data, 
respectively. In each parallelogram, we are looking down a tube into the supernova debris, 
since there is a unique correspondence, given by equation (1), between the Doppler shift of 
the observed emission and the depth in the supernova debris. To locate these parallelograms, 
we measure the differences in right ascension and declination from a given acquisition star 
and calculate the angular displacement between a pair of observations. We then shift the 
observations to a common coordinate system and measure the parallelograms directly from 
the superimposed images. For each slit orientation, we sum the pixels in the spatial direction 
(see Figs. 4 & 6) between the boundaries of the parallelogram and plot the average intensity 
(Figs. 5 & 7). We find that summing between different points, in the spatial direction, 
has a negligible effect on the intensity ratios obtained, which are smoothed (using a width 
of 5 pixels for the smoothing window) to remove anomalies due to noise or insufficient 
signal. Uncertainties in these and all subsequent measurements of reverse shock emission are 
dominated entirely by systematic errors due to contamination by emission from components 
other than the reverse shock. Although it is difficult to judge the fractional contribution of 
these contaminants, we estimate that the uncertainties are typically ~ ±10%. 
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3.1. Ha EMISSION 

We first discuss the G750L observations. Both blueshifted and redshifted emission 
features from the reverse shock are prominent in the more recent spectra on the left column 
of Fig. 4, especially just inside the inner circumstellar ring. But unfortunately, the regions 
where the reverse shock emission is brightest in these spectra usually do not overlap with the 
regions observed in the earlier spectra. For example, consider the two spectra overlapping 
in parallelogram HI. It is clear that the earlier (p8243/vl0, 1999 September 18) spectrum is 
dominated by relatively low velocity (< 3000 km s _1 ) Ha and [O i] AA6300, 6364 emission 
from the supernova core. Hence, the corresponding flux ratios are not representative of 
emission from the reverse shock. The later (pl0263/v71, 2004 July 18) spectrum is likewise 
dominated on the blue side by Ha and [O i] core emission. The very strong, blueshifted 
reverse shock emission that is evident on the north side of this spectrum does not fall in 
parallelogram HI. Redshifted emission in the velocity range 6500 < v < 8500 km s~ x on the 
southern edge of the later spectrum may be attributed to Ha emission by the reverse shock 
and also by [S ii]; this emission is not evident in the earlier spectrum. Likewise, very little 
reverse shock emission is present in parallelogram H2. 

Parallelograms H3 and H4 largely overlap and both capture the evolution of the red- 
shifted emission from the reverse shock on the southern side of the debris. In H3, the average 
intensity increases from 1999 September 18 to 2004 July 18 by a maximum factor of 2.0, 
while in H4, it increases from 1997 October 6 to 2004 July 18 by a maximum factor of 2.9. 

Parallelogram H6 captures the reverse shock only at high redshift (v > 8000 km s" 1 ) and 
high blueshift (v < —8000 km s _1 ). In these velocity ranges, the average intensity increases 
from 1997 October 6 to 2004 July 23 by factors of 1.7 and 1.5, respectively. 

Parallelogram H5 contains the only pair of Ha observations which cleanly captures the 
surface emission from the near side of the debris. The spectra shown in Fig. 5 represent 
only the emission seen in the blueshifted streaks in the respective spectra (p8243/vl0 and 
pl0263/v73). The average intensity in these spectra increased by a maximum factor of 2.7 
from 1999 September 18 to 2004 July 23. 

In summary, in the regions where we can measure the evolution of Ha emission from 
the reverse shock, both the blueshifted and redshifted average intensities have increased by 
maximum factors ~ 3, over periods of about 5 and 7 years, respectively. 
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3.2. Lya EMISSION 

The increase in average intensity is more dramatic in the G140L observations. The data, 
shown in Fig. 6, appear to be clean of contaminants and thus we believe that the intensity 
ratios plotted are reliable. 

Parallelogram LI overlays the northern side of the circumstellar ring and captures 
blueshifted emission from the reverse shock on the near side of the debris. In this region, the 
maximum average intensity has increased by a factor of 3.4 from 2002 October 29 to 2004 
July 20. 

Parallelograms L2 and L4 are partly coincident on the northern side of the debris and 
both capture much of the blueshifted reverse shock emission. In L2, we see an increase by a 
factor of 3.9 from 2001 September 24 to 2004 July 20. In L4, we see an increase by a factor 
of 9.4 from 1999 October 7 to 2004 July 20. 

Parallelograms L3, L5 and L6 give us a glimpse of the evolution of the redshifted Lya. 
The reverse shock is most prominent in L5, where we see an increase by a factor of 5.7 from 
1999 October 12 to 2004 July 20. Parallelograms L3 and L6 show increases of 9.4 and 7.5, 
from 1999 Sept 28 and Oct 13, respectively, to 2004 July 20. 

We thus see much greater increases of reverse shock emission in Lya than we do in Ha. 
For example, from September 1999 to July 2004, the blueshifted Lya in L4 increased by a 
factor of 9.4, while the blueshifted Ha in H5 increased by a factor of 2.7. In the same period, 
the redshifted Lya in L3 increased by a factor of 9.4, while the redshifted Ha in H3 increased 
by a factor of 2.0. Unfortunately, however, the regions where we can clearly measure the 
evolution of the reverse shock emission in Ha and Lya do not overlap. Therefore, we cannot 
be sure that the ratio of Lya to Ha emission from the reverse shock has increased. 

The redshifted emission from the reverse shock seen in Ha has a very different velocity 
profile from that seen in Lya (compare Figs. 4 & 6). As we discuss below, the Lya emission 
seen from the reverse shock is profoundly modified by resonant scattering, and so we cannot 
meaningfully compare the evolution of the emission in Ha and Lya from this part of the 
reverse shock. 



4. Ha VS. Lya EMISSION 

Our most recent (July 2004) STIS observations give us our first opportunity to compare 
the surface and interior emission in Ha (pl0263/v71) and Lya (pl0263/v76) through the 
same slit location and at the same epoch. Since the Ha photons are unaffected by resonant 
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scattering (M03), they are a more reliable tracer of the reverse shock intensity than the Lya 
ones. Furthermore, as discussed in §2.1, we can differentiate the surface from the interior 
emission more clearly in Ha than Lya. We construct masks to isolate the reverse shock, 
using the G750L observation as a template (Fig. 8a). We spatially sum the data and plot 
the average intensities in Fig. 9. Positions 1 to 5 are labeled on Fig. 8c; their corresponding 
positions on the G750L and G140L data are shown in Figs. 8a & 8b, respectively (assuming 
that all the emission comes from freely-streaming hydrogen atoms). 

Figure 9 shows that the ratio of Lya to Ha surface emission peaks at 13.1 on the near 
side and dips to a minimum of 2.8. From the Ha-to-Lya production ratio of 0.21 alone (M03), 
we expect 4.8 Lya photons to be emitted for each Ha photon. The fact that we observe more 
than twice this ratio is evidence of resonant scattering. Half of the Lya photons produced 
on the near side of the debris are emitted outward towards the observer. The other half are 
emitted inward towards the unshocked debris. Most of these are backscattered out of the 
debris, towards the observer. Ha photons are unaffected by resonant scattering. Hence, on 
the near side of the debris, we expect to observe a Lya/Ha ratio that is about twice the 
production ratio - i.e., a factor ~ 10. 

On the far side of the debris, the Lya-to-Ha ratio ranges from 0.71 to 11.5. We expect 
to observe a ratio less than the production ratio because most of the Lya photons produced 
there are scattered outwards, away from the observer. This is the case for Lya photons 
coming from the region labeled "1" in Fig. 8c. However, we do see Lya photons coming from 
region 2. Evidently, these photons are produced near the limb of the debris and can reach the 
observer without having to propagate through the hydrogen atoms in the supernova debris. 

5. INTERIOR EMISSION 

Besides comparing surface emission in Ha and Lya, Fig. 9 also shows the ratio of the 
surface to interior emission for both Ha (pl0263/v71) and Lya (pl0263/v76). We have 
masked out the core emission in Ha, which originates from Location 5. (Lya emission 
cannot emerge from the supernova core, owing to the high optical depth of the debris.) For 
Lya, the ratio of surface to interior emission peaks at 1.6 on the near side and 1.1 on the 
far side. For Ha, the ratios are 1.8 and 1.2 respectively. These ratios show that interior 
emission is comparable to surface emission in both cases, which means that the ratio of Lya 
to Ha interior emission tracks its corresponding surface ratio closely. Resonant scattering 
can account for the surface emission ratios discussed previously, but not for the origin of the 
interior emission (see §7). Overall, the ratio of the total Ha to Lya flux, which includes both 
surface and interior emission, is 0.17, a value consistent with the production ratio (0.21). 
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We measure the Ha surface and interior emission for observations pl0263/v72 and 
pl0263/v73 as well. These observations have the same slit orientations as pl0263/v71, but 
are positioned just west and east of it respectively (Fig. la); pl0263/71 is centered on the 
remnant. We find that, as in the case of pl0263/v71, the interior emission tracks its surface 
counterpart closely (Fig 10). On the northern (i.e., near) side of the remnant, the surface 
emission seen in the western slit is generally brighter than that in the east. But on the 
southern (i.e., far) side of the remnant, the emission seen in the eastern slit is stronger than 
that in the west. These results imply a departure of the reverse shock surface emission from 
cylindrical symmetry. We note that the X-ray (Burrows et al. 2000) and radio (Gaensler et 
al. 1997) images also show departures from cylindrical symmetry. No obvious asymmetry is 
evident in the interior emission. 



6. GLOBAL EVOLUTION OF REVERSE SHOCK EMISSION 

To construct a global light curve for the reverse shock emission in Ha, we first obtain 
profiles of the total velocity-dependent Ha flux using a combination of the 0'.'2, 0'/5 and 2" 
observations and the data reduction procedures described in §2. Again, the statistical errors 
are negligible and the main source of uncertainty comes from the systematic error, which 
stems from not knowing the shape of the emission line profile at low velocities. Michael 
et al. (2002) demonstrated, albeit in a slightly different context (i.e., in the X-ray), that 
the line profile is strongly dependent on the geometry of the system. For example, if the 
system is a cylindrical ring expanding in the equatorial plane, a "double-horned" profile 
results, with very little low- velocity emission and the "horns" centered at ±i> s , where v s is 
the freely-streaming, projected velocity of the debris in the rest frame of the reverse shock. 
By contrast, a radially-expanding spherical shell of gas produces a square-topped profile and 
a filled, uniform sphere produces a parabolic profile. 

To estimate the total line emission, we interpolate linearly over the masked data to 
obtain the Ha profiles as indicated in Fig. 11, and integrate over velocity to obtain the total 
fluxes listed in Table 2 and plotted in Fig. 12. We include the profile obtained by S05. We 
estimate very conservative lower bounds (typically ~ 45%) of the fluxes by excluding the 
interpolated profiles in the velocity ranges —3000 < v < 3000 km s _1 covered by the central 
slit. We estimate upper bounds (typically ~ 15%) by making spline fits to interpolate the 
overall profiles through the central portions. 

Since the 0'.'8 slit used in the Magellan observation by S05 did not fully capture the 
eastern and western portions of SNR 1987A, the resulting profile underestimates the amount 
of low- velocity emission actually present. Hence, it is not surprising that portions of our 
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profiles in Fig. 11 exceed the Magellan/LDSS3 one. The minimum and maximum velocities 
to which the flux extends appear to be unchanged, and they suggest that the expanding 
debris is more extended towards than away from us. We are unable to determine these to a 
greater accuracy, as the high-velocity emission is contaminated by [O i] and [S n] emission. 
Using 2" G750L data from STIS, S98 measured a total Lya flux of (1.25 ± 0.51) x 10" 12 
erg cm" 2 s" 1 on Day 3743. Since the production ratio of Ha to Lya photons is 0.21, we 
infer an approximate Ha flux of (4.86 ± 1.98) xlO" 14 erg cur 2 s" 1 . This figure is uncertain 
because the presence of resonant scattering results in an unknown conversion factor between 
the Lya photons produced and those actually observed. Nevertheless, we include this data 
point in our Ha light curve (Fig. 12), along with the (1.99 ± 0.22) xlO" 13 erg cm" 2 s" 1 result 
obtained by S05 on Day 6577. We then find that the total reverse shock flux has increased 
by a factor of about 4.1 over 2834 days (~ 8 years). 

There are two limiting cases of how the reverse shock can behave. The self-similar 

n-~ 3 

solution of Chevalier (1982) shows that the radius of the reverse shock varies as R oc t™^, 
from which we infer that the flux across the shock is 

2n — 9+4s — ns 

F(xpvR 2 (xt —s , (2) 

where the outer part of the freely-expanding supernova debris (v = R/i) and the ambient 
medium have the density profiles p oc t~ 3 v~ n and p m oc R~ s respectively. For a uniform 
medium (s = 0) and n = 9 (Eastman & Kirshner 1989), F oc t. However, the blast wave 
is now overtaking density protrusions on the equatorial ring and each encounter will send a 
reflected shock inward towards the reverse shock. If the merging of these shocks brings the 
reverse shock to a halt, then we have F oc t n ~ 4 = t 5 . Fitting a linear function, F = F Q (t — t ), 
to the data, we obtain Fq = 5.40 xlO" 17 erg cm" 2 s" 1 day" 1 and to — 2853 days. If instead 
we fit the power law, F = Fit 5 + F2, then we have F\ = 1.29 x 10" 32 erg cm" 2 s" 1 day" 5 
and F 2 = 4.87 x 10" 14 erg cm -2 s" 1 . Note that the slope of this line is determined almost 
entirely by the initial STIS observation of Lya and by the final Magellan observation. The 
intervening STIS observations of Ha are consistent with this fit, but their uncertainties are 
so great that they do not significantly constrain it, and we are unable to discern between 
the two limiting cases discussed. 



7. DISCUSSION 

Our analysis of the STIS data presented above raises more questions than it answers, 
and a detailed explanation of the observations will require theoretical modeling beyond the 
scope of the present paper. In this section, we summarize the main puzzles and suggest 
possible solutions. 
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7.1. THE Lya/Ha RATIO 

As shown in §4, resonant scattering strongly suppresses the propagation of Lya through 
the debris of SNR 1987A. Photons come into resonance when the "blueness" of their fre- 
quency is counteracted by the expansion redshift of the debris. A photon emitted in the blue 
wing of the emission profile sees an increasing optical depth and has a very high probability 
of being scattered before it penetrates far into the debris. By contrast, a photon emitted 
in the red wing becomes more redshifted as it travels through the debris, and has a good 
chance of passing through it. M03 calculated that the Sobolev optical depth for resonant 
scattering of a Lya photon is r s fa 1000. For photons emitted on the near side of the debris, 
half make the journey inward towards the supernova core. Most of these are backscattered 
towards the observer, effectively doubling the number of expected Lya from the near side. 
On the far side, the very same effect suppresses the number of Lya photons which make the 
journey through the debris and to the observer. 

M03 investigated this process using one-dimensional, "two-stream" (Schuster 1905), 
Monte Carlo simulations, but remarked that a more detailed analysis was beyond the scope 
of their study. Such an analysis will have to model the resonant scattering of Lya emission, 
in three dimensions, through an expanding debris, of characteristic thickness L ~ vt ~ 10 18 
cm, where v ~ 10, 000 km s" 1 and t — 18 years. If the path length traveled by the photon in a 
single scattering is I, then it will be redshifted by a amount u l/ct due to the expansion. The 
path length is in turn computed from the optical depth, r, generated from the probability 
for scattering, e _r . A complication arises from the fact that there is a correlation between 
I, the frequency of the photon, the temperature of the gas in the debris and the projected 
velocity of the atom it encounters (e.g., Zheng & Miralda-Escude 2002). If the Doppler 
core is involved, the surface of last scattering is located at a distance ~ 10~ 4 T^L from the 
emission surface, where T 100 = T/100 K and T is the temperature of the gas. For values 
of parameters representative of the freely-expanding debris of SNR 1987A (n fa 100 cm" 3 , 
T fa 100 K), we find that about 99% of the Lya photons are resonantly backscattered, 
resulting in a Lya-to-Ha ratio ~ 10 on the near side of the debris. The details of these 
calculations are somewhat complicated and will be presented elsewhere (Heng & McCray, in 
preparation). 

7.2. IS CHARGE TRANSFER RESPONSIBLE FOR THE "INTERIOR" 

EMISSION? 

We have shown in §5 that "interior" emission is present in both Ha and Lya and that 
it is comparable in strength to the surface emission. If we stick to the assumption that 



-13- 



the photons are emitted by freely-streaming hydrogen atoms, then we need to identify an 
emission mechanism which produces photons at a location intermediate between the reverse 
shock surface and the supernova core, with a production ratio of Lya/Ha 5. The interior 
emission, as shown in Fig. 9, is present across the broad velocity range ~ ±10, 000km s _1 . 

One possible mechanism is the interaction of the atoms with relativistic particles that 
have been accelerated at the reverse shock and have diffused upstream into the unshocked 
supernova debris. We cannot say with confidence whether this mechanism can explain the 
observed interior emission, since we have no quantitative theory to estimate the flux of 
relativistic particles in the reverse shock. 

Another possible explanation for the interior emission is that it does not actually origi- 
nate from freely-streaming hydrogen atoms in the supernova debris, but instead comes from 
hydrogen atoms resulting from charge transfer/exchange reactions of the freely-streaming 
atoms with protons in the hot gas beyond the reverse shock. The resulting atoms will have 
a broad distribution of radial velocities similar to that of the shocked protons, centered at 
~ 4000 km s -1 rather than ~ 10000 km s^ 1 . Photons emitted by subsequent excitation of 
these atoms will appear to originate from the interior of the supernova debris. However, we 
cannot be sure that this mechanism can account quantitatively for the observed emission 
without a detailed calculation, which is beyond the scope of the present paper. 

7.3. THE LIGHT CURVE OF THE REVERSE SHOCK EMISSION? 

The global evolution of the Ha emission (Fig. 12) is of great interest, since one can 
directly infer from it the net flux of atoms across the reverse shock. As S05 have pointed 
out, the ionizing radiation originating from beyond the reverse shock surface may photoionize 
the hydrogen atoms before they reach it and suppress the reverse shock emission. At present, 
the ionizing flux is increasing more rapidly than the reverse shock emission. If present trends 
continue, the reverse shock emission will vanish in < 7 years from now. 

However, as we have shown, the STIS observations do not significantly constrain this 
evolution. As described by S05, it is possible to track the global evolution of the reverse shock 
emission in Ha with ground-based observations. To do this accurately, we must obtain deep, 
low- or moderate-resolution spectra using a slit wide enough to include the entire reverse 
shock surface. This avoids the technical problems discussed in §6. 

If the Cosmic Origin Spectrograph (COS) is installed on the Hubble Space Telescope, 
it will be possible to measure both the local and global emission of Lja from the reverse 
shock. Such observations will be extremely valuable for understanding the effects of resonant 
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scattering and charge transfer as well. 

8. SUMMARY 

Using STIS data, we have tracked both the local and global evolution of the Ha and 
Lya emission from the reverse shock in SNR 1987A. The main results of our study are: 

1. The average Ha intensity has increased locally by a factor ~ 3 for both blueshifted 
and redshifted emission, over periods of about 5 and 7 years, respectively. For the average 
Lya intensity, the factors are both ~ 9 over about 5 years. However, we cannot compare 
these factors directly, as the regions where we can clearly measure the evolution of the reverse 
shock emission do not overlap. 

2. In a comparison of the emission from the reverse shock in Ha and Lya, taken through 
the same slit at the same time (July 2004), we find clear evidence that the transfer of Lya 
radiation through the supernova debris is suppressed by resonant scattering. 

3. In addition to the emission from the surface of the reverse shock, we see emission of 
comparable intensity in both Ha and Lya that appears to come from interior to the surface. 
This emission may be the result of charge transfer reactions near the reverse shock surface. 

4. In Ha observations taken with three adjacent slit positions during July 2004, we 
find departures from cylindrical symmetry in the reverse shock surface emission, with the 
brighter emission towards the northwest and the southeast. 

5. Using data from STIS, we construct the light curve of the reverse shock emission in 
Ha, which has brightened by a factor ~ 4 over about 8 years. The STIS data are consistent 
with current ground-based observations by Magellan (S05), but have much greater systematic 
uncertainties. Consequently, we are unable yet to determine the acceleration rate of the 
brightening of the reverse shock. With future ground-based observations of Ha, we should 
be able to measure this light curve with sufficient accuracy to determine whether and when 
photoionization will turn off the reverse shock emission. 
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Table 1: STIS Observations Used 



Label 


Date 


TvDe 


Slit Width 


Slit P A (dee:) 


Fxnosnre Time (si 


n1 0263/V73 


2DD4 Tnl 23 


O7 c )0T, 


0"2 


1 70 7 

1 1 J. 1 


^468 


d 10263 /v72 


2004 Jul 18 


G750L 


0"2 


179.7 


5468 


r>10263/v71 


2004 Tnl 18 

\J yj^t ' ) 1X1 X kJ 


G750L 


0"2 


1 79 7 

1 ( J. 1 


5468 


r>8243/v10 


iqqq Q Pr) f 18 

UJ J kJ LJ u X LJ 


G750L 


0"2 


-149 


0983 


n7434 /v10 


1007 Oct 6 


G750L 


0"2 


-130 


0590 


r>10263/v76 


2004 Tnl 20 


G140L 


0"2 


1 79 7 


5350 


r>9428 /v4D 


2002 Oct 20 


G140L 


0"2 


-95 3 


10600 


P 9114/v20 


2001 Sept 24 


G140L 


0"2 


-139.3 


10700 


P 8243/v23 


1999 Oct 13 


G140L 


0"2 


-124.8 


10665 


P 8243/v22 


1999 Oct 12 


G140L 


0"2 


-124.8 


10665 


P 8243/v21 


1999 Oct 7 


G140L 


0"2 


-124.8 


10478 


P 8243/v20 


1999 Sept 28 


G140L 


0"2 


-124.8 


10478 


pl0263/v74 


2004 Jul 23 


G750L 


0"5 


179.7 


5282 


pl0263/v75 


2004 Jul 23 


G750L 


0"5 


179.7 


5282 


p9428/v35 


2002 Oct 29 


G750L 


2" 


-93.8 


4832 


p9328/v2 


2002 June 8 


G750L 


2" 


110.7 


4970 


p8872/v2 


2000 Nov 3 


G750L 


2" 


-92.8 


4830 



Table 2: Total Ha Flux Obtained from STIS Observations 



Date 


Observations 


Flux 

(10~ 13 erg cm" 2 s _1 ) 


Jul 2004 
Oct 2002 
Jun 2002 
Nov 2000 


pl0263/v71+72+73+74+75 
p9428/v35 
p9328/v2 
p8872/v2 


1 Q2+ - 18 

1 yo+0.26 
L - ' °-0.78 

1 04+ - 21 
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Fig. 1.— Slit overlays for the STIS observations, (a) 0"2 G750L. (b) 0"2 G140L. (c) 075 
G750L. Only P 10263/v74 and P 10263/v75 are used in our analysis, (d) 2" G750L. 
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Fig. 2. — (a) Sample G750L observation (pl0263/v71) with the following contaminants 
indicated: The vertical stripe is due to emission from interstellar Ha and [N n] AA6548, 
6584. Vertically-aligned pairs of spots are emission from the circumstellar ring at Ha, [N il] 
AA6548, 6584, [O i] AA6300, 6364, [S il] AA6717, 6731, and He I A 6678. The horizontal 
feature located midway between the spots is Ha emission from the radioactive core of the 
supernova. A fainter such feature due to core emission in [O i] AA6300, 6364 is also evident, 
(b) Sample G140L observation (pl0263/v76) with the following contaminants indicated: N v 
AA1239, 1243 emission from a hot spot; broad damping wings from interstellar H I absorption; 
geocoronal Lya emission; interstellar line absorption due to C I AA1191, 1259 and N I A1199. 
In both cases, the insert shows the slit position overlaid on an image of the circumstellar 
ring. 
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Fig. 3. — Sample 2" observation (northern half of p9428/v35). (a) Original spectrum, (b) 
Spectrum after spatial-spectral correction is applied. The procedure uses the [O i] A6300 
line as a reference. 
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Fig. 4. — Montage of the 0'.'2 G750L observations used to track the local evolution of Ha 
emission. Each pair of observations is labeled as shown, with the newer (2004) observation 
being displayed in the left column; the corresponding slit overlays are in the right column. 
In all of the slit overlay panels, North is down. The newer (2004) observation always has 
its slit position aligned parallel to the North-South line. The region of overlap is bracketed 
by the lines shown. There are four points of overlap and summing between the innermost 
or outermost ones, in the vertical direction, has a negligible effect on the intensity ratios 
obtained (right panel of Fig. 5). 
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H1: p10263/v71 vs. p8243/v10 




v (10' km s-') 
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Fig. 5. — Left: Average Ha intensity, I V) in units of 10~ 17 erg cm -2 s -1 (km s" 1 ) -1 (arcsec)~ 2 , 
plotted vs. velocity (v). Each panel is labeled according to the scheme in Fig. 4. In all of 
the panels, the solid curve corresponds to the observation listed first. Right: Ratio of the 
new to the old intensities. 
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Fig. 6. — Montage of the 0'.'2 G140L observations used to track the local evolution of Lja 
emission. Each pair of observations is labeled as shown, with the newer (2004) observation 
being displayed in the left column; the corresponding slit overlays are in the right column. 
In all of the slit overlay panels, North is down. The newer (2004) observation always has 
its slit position aligned parallel to the North-South line. The region of overlap is bracketed 
by the lines shown. There are four points of overlap and summing between the innermost 
or outermost ones, in the vertical direction, has a negligible effect on the intensity ratios 
obtained (right panel of Fig. 7). 
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Fig. 7. — Left: Average Lya intensity, I V) in units of 10~ 15 erg cm -2 s _1 (km s^ 1 ) -1 
(arcsec)" 2 , plotted vs. velocity (v). Each panel is labeled according to the scheme in Fig. 6. 
In all of the panels, the solid curve corresponds to the observation listed first. Right: Ratio 
of the new to the old intensities. 
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Fig. 8. — (a) Ha surface emission (pl0263/v71) from the reverse shock isolated by masks, 
(b) Lja surface emission (pl0263/v76) with the same masks applied, (c) Schematic repre- 
sentation of the supernova debris with the boundary being defined by the reverse shock. For 
freely-expanding debris, there is a unique correspondence between velocity and the origin of 
the emission along the line of sight. 
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Fig. 9. — Top left: Surface (solid curve) vs. interior (dotted) average intensity, I PiV , in units 
of 10 6 photons cm -2 s _1 (km s -1 ) -1 (arcscc)~ 2 , for both Lya and Ha emission from the near 
side of the debris. Top right: Ratio of Lya-to-Ha surface emission (solid), Lya-to-Ha interior 
emission (dot-dashed), surface-to- interior emission for Lya (dotted) and Ha (dashed), for 
emission from the near side of the debris. Bottom left & right: Same as for top left and right 
respectively, but for emission from the far side of the debris. 
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Fig. 10. — Top 1st: Surface (solid curve) vs. interior (dotted) average intensity, I v , in 
units of 10~ 17 erg cm" 2 s' 1 (km s" 1 )" 1 (arcsec) -2 , plotted vs. velocity (v), for blueshifted 
emission from observation pl0263/v72. Top 2nd: Same as top 1st, but for pl0263/v73. Top 
3rd: Ratio of surface-to-interior emission for pl0263/v72 (bold solid curve) and pl0263/v73 
(bold dotted curve). Top 4th: Ratio of surface-to-surface (bold solid curve) and interior-to- 
interior (bold dotted curve) emission for pl0263/v73 vs. pl0263/v72. These ratios describe 
the east-to-west asymmetry in the emission. Bottom: Same as for the top row, but for 
redshifted emission. 
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Fig. 11. — Total Ha velocity-dependent flux, F v , in units of 10 13 erg cm 2 s 1 (km s 1 ) 1 , 
from the reverse shock. The Magellan/LDSS3 profile shown is from Smith et al. (2005). 
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Fig. 12. — Global evolution of Ha flux from the reverse shock. The thin line and bold curve 
correspond to F oc t and F oc t 5 fits, respectively. 



